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The coupling of a scalar field to electromagnetic field via the Chern-Simons term will rotate 
the polarization directions of the cosmic microwave background radiation. The rotation angle which 
relies on the distribution of the scalar field on the CMB sky is direction dependent. Such anisotropies 
will give rise to new distortions to the power spectra of CMB polarization and it can be used to 
£C) • probe the detailed physics of the scalar field. In this paper we use the updated observational data 

to constrain the anisotropic rotation angle in a model independent way. We find that the dominant 
effect of the anisotropic rotation on CMB comes from its variance and it is constrained tightly by 
O^l ■ the current data. 

USTC-ICTS-13-03 



(N 



x 



I. INTRODUCTION 



o 

As is well known the coupling of the electromagnetic field to a given external field via the Chern-Simons term 
p^AyF^ in the Lagrangian density will give rise to the rotation of linear polarization of the photons as they propagate 
in the background spacetime in which p^ does not vanish [![. Here is the electromagnetic vector field, the 
antisymmetric tensor of the electromagnetic field is — d^Ay — d^A^ and F^" = (l/2)e A " ypCT F pa is its dual. Such 
a rotation is independent of the frequencies, different from the Faraday rotation experienced by the electromagnetic 
' waves when they pass through a magnetic field. This phenomenon may be observable on large scales, for instances on 
astrophysical or cosmological scales a tiny coupling might induce a significant rotation angle potentially be observable 
by the current or near future experiments. The Chern-Simons coupling violates Lorentz and CPT symmetries in the 
background, but the gauge symmetry is preserved in cases when p^ is constant or a derivative of a scalar field, i.e., 
Pti oc dfi(f). In cosmology such a rotation will change the angular spectra of the polarization of the cosmic microwave 
background radiation (CMB), especially rotate partially the E mode polarization into the B mode one 2]. Since 
qq • both E and B modes are changed, we can use the high quality CMB observational data to detect the rotation angle 
even though no explicit B mode signal has been detected. This opened a new window to test the fundamental 
r—i . symmetries by the cosmological probe and since the first work Q it has stimulated many research interests in this 
' field [|[ . Currently the contraint on the homogeneous rotation angle from the Wilkinson Microwave Anisotropy Probe 
(WMAP) 9-year data is \ = -0.36 ± 1.24 ± 1.5 degree @. 
£C) • If the external field p M is a constant vector, there will be inconsistencies between the electromagnetic and gravita- 
tional field equations if the law of gravity is governed by Einstein's general relativity or its covariant extensions [6j. 
Alternatively the external field may be originated from derivatives of some scalar fields, e.g., p^ — (c/M)d fl <j> with M 
a certain mass scale which indicates the cutoff energy scale of the effective theory and c the dimensionless coupling 
constant. It is present in models where the scalar field <f> is a pseudo Goldstone boson produced by U(l) symmetry 
breaking. Like Axion, this scalar field has the shift symmetry and only couples to other matter derivatively. It is 
expected to interact with the electromagnetic field naturally via the Chern-Simons term. The field <f> may take the 
role of dark energy Q or not @|. It was also motivated by the leptogenesis models @ in which the baryon number 
asymmetry is produced in thermal equilibrium [To| . this violates one of the Sakharov's conditions 11 1 obeyed in tradi 



tional baryogenesis models because CPT symmetry is violated by the evolution of 4> along the time direction [l2j|. In 
this case all the equations are consistent, but the rotation angle which depends on the scalar field is not homogeneous 
across the sky because the scalar field, as a dynamical component, has different values at different positions. The 
fluctuations, treated as small effects, will bring new distortions to the CMB spectra [l3| very similar to the effects of 
weak lensing on CMB [l4[ . Measuring the direction dependence of the rotation angle with the CMB data is important 
to investigate the dynamics of the new scalar field (p. Up to now the CMB constraint on the power spectrum of the 
rotation angle from WMAP 7-year data was put in Ref.[lH based on the formalism developed in Refs. 0, H3]- They 
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obtained the result that for the multipoles from I = 1 to I = 512, there is no evidence of non-zero power spectrum for 
the rotation angle within 3a. 

In this paper we revisit this problem and compute the new constraints on the anisotropic rotation by analyzing 
the recently released WMAP 9- year data combined with data from the ground based CMB experiments QUaD and 
BICEP. We will not adopt any model for the scalar field. We find that the dominant effect on the CMB spetra of the 
anisotropic rotation angle is through the suppression factors of its variance so the current data put a strong constraint 
on the variance as well as isotropic part of the rotation angle. In comparison, current observational data are not 
sensitive to the detailed shape of the power spectrum of the rotation angle. 

This paper is organized as follows. In Sec. II, we introduce the relevant formalism for computing the rotated power 
spectra of CMB. In Sec. Ill, we present our result of the constraint on the homogeneous and anisotropic rotation 
angle from current observational data. Sec. IV is the conclusion. 



II. FORMALISM OF ROTATED POWER SPECTRA 



The full Lagrangian density of the Maxwell theory modified by the Chern-Simons term (without other sources) is 

£ = ~ F^F^+p^F^ . (1) 

This Lagrangian density is not gauge-invariant, but the action integral S = J Cd 4 x is gauge-independent because in 
this paper we only consider p^ = (c/Af)3p0 and it is easy to see that after integral by part the Chern-Simons term 
can be rewritten as — {ccj)/2M)F tlv F li,v . 

When applying the modified theory to CMB, the Chern-Simons term is expected to be significant only in the era 
after recombination. Before that, CMB photons coupled tightly to the electrons and the Chern-Simons coupling 
is negligible compared with the dominated Compton scattering. After the last scattering off the electrons, the 
polarization direction of the photon rotated by an angle [l3T ] 



X 



= -J Lss P»dx»(X) = -^A<f ) , (2) 



where the integral along the light path (parameterized by A) is from the light source to the point of observation. For 
CMB, the light source is located at the last scattering surface and one can see that the rotation angle only relies 
on the total difference of the field between the source and the observation point, A(f> = 4>(xlsStVlss) ~ 4>o, where 
xlss represents the position on the last scattering surface and r\ is the conformal time. The observation point is 
fixed but on the last scattering surface the scalar field <f> varies with the positions. In the spatially flat universe 
xlss = (va — Vlss)™, the rotation angle relies on the direction n through the dependence of the field on the position. 
Accordingly the Stokes parameters of linear polarization got a rotation 

(Q±iU) = exp(±i2x)(Q±iU) , (3) 

where we have used tilde to denote the rotated parameter. More details can be found in Ref. fl3l ]. 

On the full sky the temperature and polarization fields can be decomposed in terms of appropriate spin-weighted 
harmonic functions [l8j : 

T(n) = ^2a T ,i m Yi m {n) 

hit 

(Q±iU)(h) = £ a±2,lm ±2Yi m (n) . (4) 

ill) 

The expressions for the expansion coefficients arc 

a T ,im = j dVLY? m {h)T{h) 
a± 2 ,im = J dn ±2 Y l * m (n)(Q±iU)(n) . (5) 
Instead of ai.im and a_2,im, it is more convenient to use their linear combinations 

a-E.lra = —{a-l.lm + a-2,im)/2 

aB,lm = i(a2,lm — a-2,im)/ 2 - (6) 
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The advantage of E/B decomposition is that it is coordinate-independent and the E and B modes represent polar- 
ization patterns of opposite parity. The power spectra are defined as 

( a *X',l'm' a XAm) = X Sl'l&m'm (7) 

with the assumption of statistical isotropy. In the equation above, X' and X denote the temperature T and the 
E and B modes of the polarization field, respectively. For Gaussian theories, the statistical properties of the CMB 
temperature/polarization maps are specified fully by these six spectra. Without the rotation induced by the Chern- 
Simons coupling, Cf B = Cf B = 0. 

The Chern-Simons term has no effect on the temperature field. The rotated polarization field ([3]) can be decomposed 
in the same way except the coefficients and the spectra should be denoted by tildes. We assume the field <j> which 
induces the rotation is also a Gaussian random variable, as we have done for studies of cosmic scalar fields in other 
cases. So we also have a power spectrum for the rotation angle. As usual, we separate the rotation angle into its 
background and fluctuation, 

x(n) = x + Sx(n) , (8) 

the background part x is homogeneous across the sky and the perturbation <5%(n) as a scalar with zero mean can be 
decomposed by the spherical harmonics, 

S x (n) =J2 b imYi m (n) ■ (9) 

lm 

With assumed statistical isotropy of &z m , the angular power spectrum of the rotation angle is defined as 

{b* llm ,b lm ) = Cf5 vl 5 m , m , (10) 
which depends on the three dimensional power spectrum of the scalar field <f> at the last scattering surface, 

47TC 2 r 

°* = 1^J dlnkP ^ VLSs)jf(kA v ) , (11) 

where A77 = 770 — tjlss and ji is the spherical Bessel function. With these equations we can calculate the rotated 
spectra analytically under the expansion of the factor exp (±2i%) into Taylor series by assuming the rotation angle is 
small. Up to the second order (for spectra the fourth order), the rotated spectra of CMB was shown explicitly in Eq. 
(69) in Ref. 0. 

In this paper, instead of using Taylor expansion, we will perform the non-perturbative calculations of the rotated 
power spectra via the computations of the rotated correlation functions. This method has been used to calculate 
the lensing effect of CMB on all scales [3, [l9j], and it was shown that it is more accurate than the approach of 
Taylor expansion. The correlation function only depends on the separation of the two points and is invariant under 
displacements. In the spherical coordinate system with the origin was set to the position of the observer, we need 
only evaluate the following three rotated correlation functions for the polarization field by taking n along the z— axis 
and n' in the x — z plane at angle /3 to the z— axis 

£ + = ((Q + iU)*(n)(Q + iU)(n')) 
£- = ((Q + iU)(n)(Q + iU)(n')) 

tx = {T(n)(Q + i&){n')) . (12) 

They are related to the power spectra as 



!+(£) = ((a*E,l m -*OB,im)( a B,J'm' ^- iaB,Vm'))-iY* m {n) 2 Y Vm ,(n') 

lm,l'm' 

= +C? B hY* m (nhYt m (n') 



lm 



j2^±i(cr+cndiM, as) 



47T 
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and 



Z-W = E 2J rr^ B ~ + *c? B )*-M . (14) 



47T 

I 



= -£ ^(Q T£ + iCf*)^) , (15) 



47T 

; 



where d mk (f3) is the Wigner small d function and cos (3 = n n' . We have included the rotated TB and EB correlations. 
The reversals of the formulae (IT31) , (jT4"f and (flip)) give rise to 



c ee + £bb =2t: J ^ Z+(l3)d l 22 (p)d cos p 
d EE _ £bb = 2n [ He[f_09)]dL a2 C8)dcoBj8 



fl 



Cf " = tt / Im[e_(/3)]d'_ 22 (/3)cicos^ 



i 

Cf £ - -2tt J Reli x (p)}d l 02(P)dcosp 



C?" = -2tt y Im[a(/3)K 2 (/3)dcos/3 . (16) 

Furthermore, according to the Eq. ([3]), the rotated correlation function is related to the unrotated correlation 
function as 

e+08) = (e X p[2i( x (nO-x(n))](Q + iC/)*(Ti)(Q + i[/)(n')) 
= (exp[2i(x(n')-x(n))])£+G9) 
= cxp[-2([Sx(n')-S X (n)} 2 m + (l3) 
= exp[-4C*(0) + 4C* (/?)]£+ (/3) 

= cxp [-4C*(0) + 4C*03)] ]T -^(Cf * + Cf B )^(fi) . (17) 

/ 

where we have neglected the small effect of the correlation between the rotation angle and the unrotated polarization 
field and in the third step we used the formula that 

(e ix ) = er {x2)/2 (18) 

for a Gaussian variable a; (here it is 5\) with zero mean. The notation C x (/3) represents the two point correlation 
function of the perturbed rotation angle, i.e., 

9/ 4- 1 

C*G8) = (Sx(n)Sx(n')) = £ _^C*R(cob/?) , (19) 

i 

and C x (0) = is its variance. Substitute Eq. flTTj) into the first equation of (jTHJ), we have 

C EE + CP B = e - 4C *(°) J2 ^(Cr + C BB ) f 4' 2 (/3)4 2 (/3) e 4CX( ^dcos/3 . (20) 
Straightforwardly one may obtain that 

C EE - C BB = cos(4x) e - 4CX (°) £ 2 -^(C EE - C EB ) j\ / 22 (/3)^ 22 (/3)e- 4CX ^)deos/3 

Cf B = sin(4x)e- 4 ^(°) £ ——^—(C EE - C^) / 22 (/3)d'_ 22 (/3)e- 4CX («dcos/3 
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Cj E = cos(2x)e- 2CX W J2 2J ^ C v E f 4' 2 (/3)4o(/3)dcos/3 = Cf E cos(2 X ) e - 2CX ( ) 
Cj B = sin(2x)e- 2 <^°> £ *-±lc$ E 4' 2 (/3)4 (/3)dcos/3 = Cf * sin(2 X ) e - 2CX ( ) 



(21) 



Now we can use these formulae to compute the rotated power spectra and use the observational data together with 
the MCMC package to search for or put the constraints on \ an d Cf. 

Before that we note the spectra C E , C BB and Cf B may be decomposed into the sums, 

c l EE =c^ E +Acr , cr=c? Q B +Acr, c eb =c eb +ac eb , ( 22 ) 

where 

C EE = cos 2 (2x) + C BB sin 2 (2x)] e - 4CX (°) 

C BB = [Cf E sin 2 (2x) + G BB cos 2 (2x)] e - 4CX (°) 

= i sin(4x)(Cf E - C BB )e- iCX{0) (23) 

together with CT E and C^ 3 they depend on the power spectrum of the rotation angle Cf only through the variance 
C x (0). Other parts rely on the correlation function of nonzero angular separation, 



AC EE + AC BB = 2tt J A£ + (/3)4 2 (/3)dcos/3 
AC EE - AC BB = 2ir J Re[A|_(/3)]4 2 (^)dcos^ 



ACf B : ,t j M*Z-<M4nW<x»P . ( 24 ) 



where 



A|_(/3) = e 4 ^- cx (°))(e- 4 ^(« - 1) £ ^(Q ££ - Q BS )^ 22 (/3) . (25) 



III. RESULTS 



In the above section, we have calculated the rotated CMB power spectra based on the full sky formalism. The results 
showed that besides the effects of homogeneous rotation x, the anisotropic rotation angle brought new distortions to 
the CMB spectra, indicated by the variance C x (0) and two point correlation function C x ((3) in the formulae ([20| and 
(1211) , both of which depend on the angular power spectrum Cf and in principle can be constrained by the data from 
CMB polarization experiments. In order to constrain the homogeneous and anisotropic rotation angle in a model 
independent way, that means we will not assume any model for the scalar field (f>, we first modify the publicly available 
MCMC package CosmoMC [20] to include the distortions in Eqs. ([2"U]) and (|2"T) by treating x as a free parameter and 
Cf as a positive otherwise undetermined function of I, and then perform a global fit to the CMB data. This method is 
different from that developed in [15l |16L |21| , which adopted a quadratic estimator to constrain the anisotropic rotation 
angle. The main difference is that in our method we treat the rotation angle as a Gaussian random field, its statistics 
is isotropic, so that the rotated polarizations of CMB which were identified as the observables are also statistically 
isotropic and there will be no off-diagonal correlations between different multipoles I and I'. However, the quadratic 
estimator method considers the ensemble of CMB rotated by a fixed rotation angle x(^) which is not averaged and 
breaks the statistical isotropy of the rotated polarization field of CMB. Hence the off-diagonal modes of the CMB 
correlations do not vanish and depend explicitly on the rotation angle x(^)j since the averaging is done only on the 
CMB ensemble. Especially, in the quadratic estimator method, some of the CMB correlations contain 0(Sx) terms, 
based on which the quadratic estimator is constructed. Different from this, in our method both the unrotated CMB 
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field and the rotation field are considered as Gaussian random fields and have isotropic statistics. The averaging are 
done on both the CMB and rotation field ensembles, only the diagonal terms are nonzero and the results do not 
depend on $x in the coordinate space, but rather on its auto-correlation functions (or its power spectrum Cf). This 
point can be understood in detail by comparing the CMB EB correlation appearing in [la ] (equation (9)) and that in 
our paper. The CMB EB correlation in [16| has off-diagonal modes, and its first term is O(x), while the second term 
is 0(8x)- In our paper, the CMB EB correlation only has diagonal Cf B modes and only 0(x) and 0(Cf) terms 
appear in Cf B (the 0(x) term is the same as that in |16j if x is replaced by :J^=). So in our analysis it is necessary 

to model the shape of Cf, in the quadratic estimator method it does not need to do so. 

In our computations, besides the WMAP 9yr temperature and polarization data [22[ (including the \ow-l TB, low-Z 
EB and high-Z TB data), we also adopt the QUaD data [23| and the BICEP data (24j. As for the cosmological model, we 
use a flat ACDM Universe with a cosmological constant. The free cosmological parameters are u> c , ©s> t, n 8 ,r, A s ), 
where w& = f^/i 2 and u> c = VL c h 2 are the baryon density and cold dark matter density respectively, O s is the ratio 
of the sound horizon to the angular diameter distance at decoupling, r is the optical depth to reionization, n s is 
scalar spectral index, r is the tensor to scalar ratio of the primordial spectrum and A s defines the amplitude of the 
primordial scalar spectrum. 

Compared with the choice for the cosmological parameters, it is difficult to choose the parameters for the anisotropic 
rotation angle. From Eqs. (j2"UI (2TJ [231 (Ml) we note that the distortions to the CMB power spectra brought by 
the anisotropies of the rotation angle can be divided into two parts. The first part appears in C EE ,Cf§ \Cf^ ', 
Cf E and Cj B , in which the anisotropies of the rotation angle suppress the CMB angular power spectra through the 
variance C x (0) via the factors exp(-4C x (0)) and exp(-2C x (0)). The second one is ACf E , AC BB and ACf B , which 
depends on the detailed shape of Cf. Generally speaking, the first one dominates the CMB power spectra distortions. 
Therefore, the data will give a strong constraint on C x (0) but a weak constraint on the shape of Cf. Furthermore, the 
constraint on C x (0) = J2i TF^* coming directly from the Cf E , Cf^ ', Cf B , Cf E and Cj B terms is larger than that 
coming directly from the AC EE , AC BB and ACf s terms. Due to this, it is difficult to obtain a precise constraint 
on the shape of Cf. We can only get a loose constraint on it. The method we adopted is to let C x (0) runs freely with 
Cf . More specifically, we adopt six binned Cf parameters Cf x , Cf 2 , Cf 3 , Cf A , Cf 5 , Cf 6 , which are the average values 
of Cf in the multipole regions [2, 100], [101,200], [201,300], [301,400], [401,500], [501, l max ] respectively (l max = 2550 
is the largest multipole of Cf we used), to model the shape of Cf. So in this case, the general parameter space is 

P = {uo bl u c , 9 S , r, n s ,r, A s ,x, C x (0), Cf l7 Cf 2 ,Cf 3 ,Cf A , Cf 5 ,Cf 6 }. (26) 

Besides these quantities we also calculate six derived parameters Cf(0) = X^=i TF^fi where i = 1 — 6, which 
can be used to compare C x (0) with Y^i=i Cf (0) ■ When running the program, we set a physical prior on C x (0) 
and Cf i7 that is C x (0) > and > 0, but we do not require C x (0) > £i=i Cf(0). The constraints on these 
parameters are illuminated in Figure [TJ from which no positive signal of nonzero CMB polarization rotation or 
possible anisotropies are found, but a strong constraint is obtained. More specifically, at 2a level of confidence, the 
constraints are —0.043 < x < 0.034 (throughout this section, radian is adopted as the unit of the angle), C x (0) < 0.020, 
Cf A < 1.1 x 10~ 5 , Cf 2 < 2.8 x 10~ 6 , Cf 3 < 2.0 x 10~ 6 , Cf 4 < 1.8 x 10" 6 , Cf b < 1.5 x 10~ 6 , Cf 6 < 5.2 x 10~ 8 , 
Cf(0) < 0.0090, Cf (0) < 0.0068, Cf (0) < 0.0081, Cf(0) < 0.0099, Cf{0) < 0.011 and Cf{0) < 0.034 respectively. 
Compared with the results presented in [l5| . it seems that the constraint obtained here is tighter. However, it is 
difficult to directly compare them with each other, since the two methods differ in several aspects. First, as we 
mentioned at the beginning of this section, we treated the anisotropic rotation angle as a Gaussian random field with 
isotropic statistics and use the rotated power spectra to constrain its anisotropies, while in [l5| the authors adopted 
the quadratic maps to estimate the fixed anisotropic rotation angle. Second, we set a positive prior on the C x (0) and 
Cf. Third, the bin size we adopted is Al ~ 100, which is different from Al ~ 50 used in [15j . 

At the same time, one can see that the constraint on C x (0) is much stronger than Yl%=i C?(0); which indicates that 
the distortion of the CMB power spectra induced by the anisotropies of the rotation angle is indeed dominated by 
the variance C x (0) which appears as suppression factors in C EE , C BB , C EB , Cj E and Cf B . Due to this, we consider 

another simpler way, in which the effects of AC EE , AC BB and AC EB are ignored and the parameter space becomes 

P = {u b , u e , Q s ,t, n„r, A s ,x, C x (0)}. (27) 

The results are displayed in Figure [5] and the 2er constraints are —0.071 < x < 0.057, C x (0) < 0.014. The constraints 
are similar to that obtained with the first way, but the constraint on x becomes weaker and the constraint on C x (0) 
becomes stronger. 

Based on the above analysis, we obtain the following conclusions about the results: (1) Currently the data put 
strong constraints on isotropic rotation angle x as weu as the variance of the anisotropic rotation angle C x (0). We 
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FIG. 1: One-dimensional constraints on the isotropic rotation angle x> the variance of the anisotropics of the rotation angle 
C x (0), binned Cf parameters Cf. and derived parameters Cf(G). In obtaining this result, we let C x (0) runs freely with the 
binned Cf parameters. 

can conclude conservatively that C x (0) < 0.020; (2) Although a precise constraint on the shape of Cf is not obtained, 
one can still give some weak constraints on it. 
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FIG. 2: One-dimensional constraints on the isotropic rotation angle \ an d the variance of the anisotropies of the rotation angle 
C x (0). In obtaining this result, we ignored the effects of ACf E , ACf B and ACf B , which depend on the detailed shape of 



IV. CONCLUSIONS 

Detecting the rotation of the CMB polarization induced by the Chern-Simons coupling is now considered as an 
important way to test Lorentz and CPT symmetries in the community. Usually people only focus on the isotropic 
rotation. Theoretically if the rotation is induced by a fixed vector, when including the gravity, the equations of motions 
are not consistent. Alternatively if it is induced by the scalar field coupling, the rotation angle must be direction 
dependent and we should include the anisotropy of the rotation angle in the more consistent analysis. Furthermore, 
measuring the anisotropy of the rotation angle is equivalent to determine the distribution of the scalar field on the 
last scattering surface and this is very important for us to probe the detailed dynamics of the scalar field. In this 
paper, in a model independent way, we studied the constraints on the anisotropic rotation angle by analyzing the 
recently released WMAP 9-year data and the data from QUaD and BICEP experiments. We found that even though 
current experiments have little to say about the shape of the power spectrum of the rotation angle, but they put 
strong constraint on its variance and its mean value (the homogeneous part of the rotation angle). Our result showed 
no significant evidence for a nonzero rotation angle. 
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